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CSj ! Abstract 

> 

Q\ . The number of BPS bound states of D-branes on a Calabi-Yau manifold depends on two 

sets of data, the BPS charges and the stability conditions. For DO and D2-branes bound 
to a single D6-brane wrapping a Calabi-Yau 3-fold X, both are naturally related to the 
Kahler moduli space A4(X). We construct unitary one- matrix models which count such 
BPS states for a class of toric Calabi-Yau manifolds at infinite 't Hooft coupling. The 
matrix model for the BPS counting on X turns out to give the topological string partition 
function for another Calabi-Yau manifold Y, whose Kahler moduli space A4(Y) contains 
two copies of Ai(X), one related to the BPS charges and another to the stability con- 
ditions. The two sets of data are unified in Ai(Y). The matrix models have a number 
of other interesting features. They compute spectral curves and mirror maps relevant to 
the remodeling conjecture. For finite 't Hooft coupling they give rise to yet more general 
geometry Y containing Y. 



1 On leave from University of Amsterdam and Soltan Institute for Nuclear Studies, Poland. 



1 Introduction 

The topological string theory has deep connections to a variety of BPS counting prob- 
lems in string theory [HE]. In this paper, we focus on the generalized Donaldson-Thomas 
(DT) invariants, namely the numbers of DO and D2 bound states on a single D6 brane 
wrapping a Calabi-Yau 3-fold X. The DT invariants are background dependent. As we 
vary the Kahler moduli of X and cross a wall of marginal stability, the numbers can jump. 
To count BPS bound states, we have to specify the stability condition, i.e. the chamber in 
the moduli space where we perform the counting. Thus, the DT invariant depends on two 
sets of data, the BPS charges and the stability conditions. In particular, the commutative 
DT invariants are defined in the chamber corresponding to the infinity in the Kahler mod- 
uli space, while the non-commutative DT invariants are defined in the chamber containing 
the origin. 

It is convenient to introduce the generating function -Zbps of the DT invariants Q a ^(n), 

Z B Ps(q,Q;n) = J2 n -An)q a Q P , (1-1) 

a,/3 

where a G Z is the DO brane charge, /3 G H2(X, Z) are the D2 brane charges, and n is a 
set of parameters which specify the chamber in the Kahler moduli space. In this paper we 
consider toric Calabi-Yau manifolds without compact 4-cycles, see figure For a manifold 
X in this class, it was shown in [3j that -Zpps is given by a certain reduction of the square 
of the topological string partition function Z top (q, Q), 

ZBPs(q, Q] n) = Z top (q, Q) ■ Z top (q, Q~ l ) . (1.2) 

reduction at n 

In this case, Z t0 p{q, Q) is expressed as a product in the harmonic oscillator form. The 
reduction means dropping an appropriate set of harmonic oscillator factors from |Z t0 p| 2 
corresponding to D0/D2 states that do not bind with the single D6 brane in the chamber 
n. 

Both Q and n are related to the Kahler moduli space M.(X) of X. The relation of n to 
the moduli space is clear since it specifies a chamber in M.(X). It is also natural to identify 
Q = e~ l in (II. ip with t being flat coordinates of M.(X) since the BPS charges couple to 
the areas of the corresponding homology cycles, i.e. the Kahler moduli. However, these 
two data appear asymmetrically in ( II. 2p . In this paper, we will present another connection 
of -Zbps to the topological string theory, in which they are treated more symmetrically. 
We will show that there is another Calabi-Yau manifold Y, whose Kahler moduli space 
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M.(Y) contains two copies of M.(X), and the topological string partition function for Y is 
related to -Zbps f° r X. For example, when X is the resolved conifold with dime M. (X) = 1, 
the corresponding Y is the suspended pinch point (SPP) geometry with dime Ai(Y) = 2. 
Similarly, when X is C 3 /Z,2, the corresponding Y is C 3 /Z^. 

We will find this relation by constructing the unitary one-matrix model whose partition 
function Z matrix (g, Q; n) is related to Zbps(?, Q', n). In particular, Z matrix is equal to Zbps m 
the non-commutative chamber (n = 0) and is equal to Z top (X) in the commutative chamber 
(n = oo). To derive the matrix model, we start with the crystal melting model [U EI] to 
count the generalized DT invariants, and use the vertex operator formalism [6, 7], in which 
the partition function is expressed as correlators of exponentials of fermion bilinears. The 
correlators are defined for all chambers in the Kahler moduli space, and we can transform 
the computation into unitary matrix integrals. This construction is closely connected to 
the free fermion picture for the topological string and Seiberg-Witten theory developed 
in [HI El ITU] . Equivalently, we can also express the partition function as a sum over non- 
intersecting paths following and generalizing [TTJ, which gives yet another derivation of 
such matrix models. 

One interesting feature of our matrix model for the conifold, in the commutative cham- 
ber, is its close relation to the so-called Chern-Simons matrix model of [T2"| IT3] . In the 
commutative chamber in our model, Q is the only parameter, and it appears only in the 
potential. The Chern-Simons matrix model also depends on a single parameter, which 
is the 't Hooft coupling. It turns out that these two parameters play the same role in 
the partition function in both models. Moreover one can consider a model with non-zero 
values of both these parameters. From this viewpoint, departure from the commutative 
to arbitrary chamber can be interpreted as turning on yet another parameter. In general 
one can consider simultaneously non-zero values of all three parameters: Q, chamber de- 
pendence, and 't Hooft coupling. This gives rise to the spectral curve encoding yet more 
general Calabi-Yau manifold Y, which contains the manifold Y described above. When X 
is the conifold, Y is a symmetric resolution of C 3 /Z2 x Z2, while Y is the SPP geometry 
as we mentioned in the above. Such Y can in principle be constructed for any initial toric 
manifold X. 

As a bonus of our matrix model construction, it sheds new light on the remodeling 
conjecture. It has been conjectured in [H] that the topological string partition function 
for this class of Calabi-Yau manifolds is completely characterized by the recursion relations 
of [T3], applied to the curve which should be identified with the mirror curve of a given 



manifold. Such recursion relations would arise if we had a matrix model formulation of 
the topological strings. In this paper we provide a construction of such matrix models in 
several instructive cases, and verify that to the leading order their spectral curves agree 
with relevant mirror curves, which is an important step towards a proof of the remodeling 
conjecture. We expect that application of our methods should lead to analogous results in 
the general case of toric manifold without compact 4-cycles. 

We also note that, for the case of C 3 , a similar approach was presented in |TT]. For 
an earlier related work, see [16]. Matrix models for other Calabi-Yau manifolds in the 
commutative chamber were derived from the topological vertex formalism or Nekrasov 
partition functions in [T7| [TSj [TjJJ |2DJ. In the course of this work we received the paper |21] . 
in which matrix models are derived in the commutative chamber also from the topological 
vertex perspective. Related ideas have been considered in [22| 123]. 

This paper is organized as follows. In section 2 we introduce matrix models for BPS 
counting and explain how they are related to the DT invariants. In section 3 we examine 
spectral curves of the matrix models and identify the corresponding Calabi-Yau geometries. 
In particular, when X is the resolved conifold, we also identify the total geometry Y for 
finite 't Hooft coupling, and discuss its relations to the Chern-Simons matrix model. The 
derivation of the matrix model is given in section 4. We end with summary and discussion 
on future research directions in section 5. 

2 Matrix Models 

In this section, we will present matrix models which count the DT invariants, namely 
the number of BPS states of DO and D2-branes bound to a single D6 wrapping a Calabi- 
Yau manifold X. In general these are matrix models for unitary matrices of infinite size, 
and arise from crystal melting interpretation of BPS generating functions. The derivation 
of these matrix models will be given in section 4. 

2.1 C 3 

When X = C 3 , the generating function of BPS invariants is given by the MacMahon 
function which counts plane partitions. We find that this BPS generating function is equal 
to the partition function of the matrix model given by 



^matrix I 



(q) = / dU det e(U\q), (2.1] 



where the integral is over the unitary group U(N), and we are interested in the limit of 
N — > oo. The integrand is given by the theta-product, 

oo 

Q(u\q) = J](l + uq k ){\ + «-y +1 ). (2.2) 

fc=0 

To perform the integral (12. ip . it is convenient to diagonalize U = diag(iti, ...,u n ) and to 
consider the integral over eigenvalues Ui = e 1 ^ as 

^matnx(g) = / ]J d<f>i 6(e*|g) JJ(e** - e^)(e-^ - e~^). (2.3) 

i i<j 

As usual, the two factors of the Vandermonde determinant come from the integral over 
off-diagonal elements of U. To perform the integral (12. 3p over eigenvalues, we expand the 
integrand in powers of q, 

9(e^|g) = 1 + e* + (1 + e"^ + e^ + e 2 ^) g + (2 + tT* + 2e** + e 2 ^) g 2 + • • • , 

and pick up appropriate combinations of e ±4< ^ 's from the measure factor in (12. 3p to cancel 
the (^-dependence in Q{e l ^\q). In this way, we can directly verify that the integral gives 
the MacMahon function, 

oo 1 

Z ma trix(g) = 1 + q + 3q 2 + Qq 3 + 13g 4 + = J[ (2.4) 

fe=i ^ ^ ' 

This is indeed the generating function of plane partitions and reproduces the counting of 

the DT invariants on C 3 if we identify the power of q as the DO brane charge. In this case, 

there is no distinction between commutative and non-commutative chambers. 

To relate this to the Chern-Simons matrix model, we make the identification of q = e~ 9a , 

where g s is the string coupling constant. For small g s , the modular transformation of 

with respect to g s gives 



,2 , 1 



Q(e l(p \e- 9s ) = e"5ii • 1 + 0(e^^) . (2.5) 



If we ignore non-perturbative terms in g s , this is equal to the integrand for the unitary 
Gaussian matrix model derived from the Chern-Simons theory on the conifold |13| . In 
fact, ( 12. ip itself has also been proposed for the topological string theory on the conifold in 
[16| . whose approach is a special case of our fermionic derivation applied to C 3 as we will 
see below. The Kahler moduli T of the resolved conifold is given by the 't Hooft coupling, 

T = g s N. (2.6) 



We are interested in the iV — y oo limit for fixed g s , namely T — y oo. It is shown in 
that the model ( 12. 3p with finite A" has an interpretation of counting plane partitions 
in a container with a wall at position N. As we will discuss in the next section, finite 
't Hooft parameter has similar wall interpretation in our more general models. From 
this perspective, A" — y oo limit in C 3 model corresponds to computing all plane partitions. 
This limit suppresses instanton corrections on the conifold, leaving only contributions from 
constant maps. For a general Calabi-Yau manifold, the sum over constant maps gives the 
MacMahon function to the power of x/2, where \ is the Euler characteristics of the Calabi- 
Yau manifold. Since x — 2 for the resolved conifold, we find that the A" — y oo limit gives 
one power of the MacMahon function, reproducing ( 12. 4p . 

2.2 Conifold 

The Kahler moduli space of the resolved conifold is complex 1-dimensional, and it is 
divided into chambers parametrized by an integer n, which is the integer part of the B- 
field flux through the P 1 [3 J . The non-commutative chamber corresponds to n = and the 
commutative chamber is at n = oo. 

We find the following matrix model in the non- commutative chamber: 

^matrixte, Q;n = 0) = JdU det ( J^jrV ) , (2-7) 

where 

Q = e-* (2.8) 

keeps track of the D2 brane charge. By expanding the integrand in powers of q and by 
performing the integral over U(N) in the N —y oo limit as in the previous example, we can 
verify that 

Z matI Uq, Q\ n = 0) = 1 + (2 - Q- 1 - Q)q + (8 - AQ' 1 - AQ)q 2 + ■■■ 

fr (l-Q<ffll-Q-W (29) 



fc=l 



'l- q k)2k 



This reproduces Zbps{(1, Q'i n — 0) in the non- commutative chamber. 
For a general chamber, the BPS partition function is given by 



The free fermion expression for Z^ps(q,Q;n), discussed in section 4, gives rise to the 
following matrix integral, 

Z mat rUl,Q;n) = f dU det ( q(qu\L ft 1 + Q~ X U- X q k )\ ■ (2-11) 

The BPS partition function and the matrix model partition function are related as 

Z B ps(<?, Q\ n) = C n - Z matrix (q, Q; n), (2.12) 

where the prefactor C n is given by 

1 -^ fl-Q- l q k 



^n^n w ■ < 2 - 13 ' 



q , 

fe=l v ^ ' k=n+l 

We also verfied (j2.12p by expanding matrix model integrand and integrating it term by 
term. The origin of the prefactor C n will be explained in section 4. Note that this prefactor 
is trivial in the non- commutative chamber, C n= o = 1. 

It is known that the BPS partition function in the commutative chamber and the 
topological string partition function are identical, up to one power of the MacMahon 
function, 

oo 1 

^BPs(g, Q; n = oo) = Z top (q = e~ 9 % Q = e~ l ) ■ J] (2.14) 

fc=i ^ ^ ' 

Since the prefactor C n reduces to the MacMahon function in the commutative limit, 



C. 



oo 1 



q 
fc=i v H 



the matrix model partition function gives precisely the topological string partition function 
in the commutative chamber, 

*„ x.. f Jr r ., * /tt C 1 + *W "+ U~Y +1 



J \k=0 V 



mm.-!* /"''III ,1 + QUqk) 

= Z top (q,Q). (2.16) 

In this way, the matrix model partition function Z matrix (q, Q; n) interpolates between Zbps 
in the non- commutative chamber and Z top in the commutative chamber. 



2.3 C 3 /^2 

Another toric Calabi-Yau manifold with dime Ai(X) = 1 is C 3 /Z,2. The matrix model 
for the non-commutative chamber is given by 



J matrix I, 



(q, Q-n = 0)= dU det (e(U\q)Q(QU\q)) . (2.17) 



For a general chamber, we can write the explicit product form of the BPS generating 
function as a matrix integral 



Z B Ps(q,Q;n) = n(l-g*)-*(i-Q?*)- fc (l-Q-V 



I— 1 „n+k\— n— k 
x — y j ^x — vy ) K*- — ^ 

fe=l 

n 'J vnLi(i+Q-^-v) ' ( 8) 



with the prefactor 

Cn = n a _ „k) k n (fi- a ^fi-Q-v)j • (2 - i9) 



[(i-i k ) k k i n \ 1 W-Q k )(i-Q- 1 <i k 

This can be verified explicitly by expanding both sides of (I2.18P in powers of q. 

Again, in this case, we have C n= o = 1 and C n=00 = Yikfi ~ Q ) ■ Thus, the matrix 
model partition function interpolates between Zbps i n the non- commutative chamber and 
Z top in the commutative chamber, 

Z mat ri x (g, Q\ n = 0) = Z top (q, Q) ■ Z top (q, Q' 1 ) = Z BPS (q, Q; n = 0), 
Zmatrix(<?, Q] n = oo) = Z top (q, Q). 

2.4 General Toric Calabi-Yau Manifold 

A toric Calabi-Yau 3-fold X without compact 4-cycle consists of a chain of P^s, which 
is resolved either by 0{— 1, —1) or 0{— 2, 0). The topological string partition function for 
such a Calabi-Yau manifold is given by flQ |2"5] 

(oo \ X/2 oo 

iitt^^ n n^-a-'-Qirt"""'*, (2.21) 
fc=i ^ y ' / i<j<i<x-ifc=i 

where x is the Euler characteristics of X, the number of P^s is (% — 1), and Q 1; ..., Q x -i are 
the Kahler moduli that measure their sizes. Depending on whether the i-th P 1 is resolved 
by C(-l, -1) or 0(-2, 0), we set s, = -1 or +1. 



The BPS partition function in the non-commutative chamber is given by 

Z B ps(q, Q; n = 0) = Z top (q, Q) ■ Z top (q, Q' 1 ). (2.22) 

This is reproduced by the matrix model partition function, 

r x_1 

^matnx(g, Q\ n = 0) = / dU det n e(siQi - • - ^Qit/lg) 31 "^. (2.23) 

■^ j=l 

Following the procedure described in section 4, it is possible to write down matrix models 
for other chambers. However, we have not attempted to derive a closed-form expression of 
the matrix model potential for a general chamber. 

3 Spectral Curves and Geometric Unification 

The eigenvalue distribution of the large N matrix model is controlled by a spectral 
curve. In particular, the resolvent is a one-form on the curve and the large N effective 
action is evaluated by its period integral on the curve. It has been argued from several 
viewpoints that spectral curves of matrix models arising from the topological string theory 
should be related to the geometry of the corresponding Calabi-Yau manifold |14[ [TS| |2"E] . 

In this section, we will identify the spectral curves and the corresponding Calabi-Yau 
geometries Y for the matrix models defined in the previous section. These geometries arise 
in the limit of infinite 't Hooft coupling. In a nontrivial case of X ^ C 3 , they contain 
two copies of the initial Calabi-Yau manifold X for a generic chamber. For the conifold 
case, we will analyze in detail yet more general geometry Y which arises for finite 't Hooft 
coupling, as well as reveal close relation between conifold matrix model in the commutative 
chamber and so-called Chern-Simons matrix model [T2| 

3.1 C 3 

As a warm-up exercise, let us describe the unitary Gaussian model, discussed in 

HS1E7], 

Z m ,tnM = f II d ^ e ~^ ? 11^ - e ^)(e-^ - e~^). (3.1) 

i i<j 

Since c/>;'s are periodic variables, it may appear unnatural to have the non-periodic poten- 
tial, e 2 3s v i. In our construction, it is the g s — > limit of the periodic integrand given in 



(12. 3p . The integrand Q(e l ^\q) has a series of zeros at = ikg s with k e Z, which becomes 
a branch cut along the imaginary axis in the limit g s — > 0. 

The spectral curve for the unitary Gaussian matrix model is given by the equation 

P3JEI] 

e x + e y + e x - y ~ T + 1=0, (3.2) 

where T = Ng s is the 't Hooft coupling. The corresponding Calabi-Yau manifold is the 
mirror of the resolved conifold. In the limit of T — > oo, the curve reduces to 

e x + e y + 1 = 0, (3.3) 

which is the mirror of C 3 . This result also arises as a special case Q,e~ T ,fi — > of a 
derivation of the conifold curve presented in the next section. 

3.2 Conifold 

In this section we analyze conifold matrix model in all chambers. From the form of 
the spectral curve and for finite 't Hooft coupling we identify the total manifold Y to be a 
resolution of the C 3 /Z2 x Z2 orbifold. In the limit of infinite 't Hooft coupling, Y reduces 
to the susspended pinch point (SPP) geometry Y, which contains two copies of the initial 
conifold geometry. 

To examine the g s — > limit of the matrix model for the conifold, let us look at the 
integrand of ( 12. lip . We find it convenient to choose the freedom of renaming U — > U~ l 
described in section 14.1.21 and consider an equivalent integrand 



m-'k) A n , n - 1Unk] = n (i + u-Y){i + ug k+1 ) 

e{QU-i\q) l£ "" V q ] 11 (l + QU-iq k )(l + Qr-We-rqW) 



(3.4) 
^lit^ -q~){L-r*4 ~ue q~-) 

Here and in what follows we set 



t = ng s . 

In order to retain interesting dependence on the chamber parameter n, we should take the 

limit g s — > in such a way that r is held finite. By using the identity, 

00 1 

log ]J(1 + Uq k ) U 2 (-U), (\g a \ « 1), (3.5) 

fc=i 9s 

where Li2 is the dilogarithm function, the integrand (13. 4p can be approximated by e~^ *■ ' 
with 

V{U) = TlogU + U 2 (-U) + Lizi-U- 1 ) - U 2 {-QU- 1 ) - L^-Q^e^U), (3.6) 



where we also took into account the shift ( 1A.1|) of the potential which arises from the 



transformation of the measure to the form which includes the Vandermonde determinant. 

Therefore 

.__ T - logg + Q) + log (1 + J*,) 

d v V = . (3.7) 

Let us define the resolvent u(u) by 

where the expectation value is taken over the large N eigenvalue distribution of U. As 
we expect to find a genus curve, we postulate the existence of a one-cut solution. In 
this case, in the weakly coupled phase of a unitary matrix model, the resolvent can be 
computed using the standard Migdal integral 

M() = ±I±d v V(v)y/ (u-a + )(u-^_ 
v; 2TJ 2m u-v ^{v - a+ )(v - o_) 

where the integration contour encircles counter-clockwise the endpoints of the cut a±. We 
perform this computation in appendix A and find 

/x J_, ( \/{cl+ + Q) (a- ~ jj) ~ y/(a- + Q) (a+ ~ u ) u + Qe T e T l 2 
A { " ' ~~ uT ° g ^ y /( a+ + Qe T )(a--u) - ^{a_ + Qe T )(a+ - u) u + Q Q^e^ 2 

(3.10) 

This form of the resolvent already takes into account the boundary condition 

u(u — y oo) ~ — . (3-11) 

u 

This condition also gives rise to two equations on the location of a± 






\/(a+ + Q)a_ - y/(a- + Q)a + ^ _ _ (r+T )/2 



a/(o + + Qe T )a_ - y/(al + Qe T )a + 
With some effort these equations can be solved in the exact form 

2 (1 - //)(! - M e 2 ) + (1 - Q)(l + /ie 2 - 2/0 



Qt e -l^'J/2. (3.13) 



o ± = -1 + e' 



(l-/ie 2 ) 2 



, „. y / (l-g)(l-e 2 )(l-^)(l-Q/i6 2 ) 
' C (l-//e 2 ) 2 
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(3.14) 



see figure [TJ The parameters fi and e are related to the chamber number r and the 't Hooft 
parameter T as 

/i = g- 1 e - r , e = e~ T/2 . (3.15) 

From the form of (13. 14j) it is clear that in the saddle point approximation the cuts are 
deformed and do not lay on a unit circle, but (as often happens in similar situations) on 
arcs which are deformations thereof. One can also verify that u(u) given by (I3.10p is indeed 
a solution to the Riemann-Hilbert problem, 

ldV 



lo+(u) + UJ-(u) 



T du 



(3.16) 



where u± are the values of u(u) right above and below the branch cut. From the resolvent 
one can also find the eigenvalue density (for fi < 1) 



p{u) 



U) + (U) — u~{u) 



J_ / (1 + fie 2 )u + 1 + Qe 2 - (1 - fie 2 ) y/(u - a+)(u -ajj 
u? ° g l (1 + ^)u + 1 + Qe 2 + (1 - /ie 2 ) v^T^ a+){u - a_) 




Figure 1: Behavior of cut end-points a + (in blue) and a_ (m redj g'zwen in 
^3. 14 ), for fixed e,Q and varying \x. For \x < 1, end-points a± are complex 



conjugate to each other. For \x = 1 we have a^ 



-1 - 



1-e 2 



and the 



cut shrinks to zero size. For \x > 1 both a± are real and spread in opposite 
directions. 



To identify the spectral curve we note first that the non-trivial part of the resolvent 
takes the form (for \l < 1) 

N (3-17) 



U)(U) 



^ l0g 



u 



l + Qe 2 1 



fie 



l + fie 2 1 + fie 2 
11 



\/(u- a + ){u - a_)). 



After identification x = uTu)(u), and setting u = e y , we find that e x and e y satisfy a 
polynomial equation. Appropriate constant shifts of x and y transform this equation into 
the following form 

e x+y + e x + e y + Q 1 e 2x + Q 2 e 2y + Q 3 = 0, (3.18) 

where 

2 1 + fJ.Q 



Qi = e 



;i + ^e 2 )(l + Qe 



■> 



1 + Qe 2 . 

(1 + nQ)(l + /ie 2 ) 



Q 3 = Q 



1 + fte 



2 



;i + e 2Q)(i + / iQ)- 

The above equation represents the spectral curve we have been after. It is interesting that 
the curve (I3.18P is symmetric under exchanges of Q, fl — Q~ x q n and e 2 = e~ T . Namely, 
the original Kahler moduli Q of the resolved conifold, the chamber parameter n and the 
't Hooft parameter T appear symmetrically in the spectral curve. We also note that the 
above form of the curve, as well as the density and the resolvent given in (|3.17j) . are valid 
for \fi\ < 1. For |/x| > 1, an appropriate analytic continuation is required. 

The corresponding Calabi-Yau manifold Y is a resolution of the orbifold C 3 /Z2 x Z2. 
There are two such resolutions, the symmetric one (also known as closed topological vertex) 
and the asymmetric one. Both of these resolutions consist of three P^s and are related to 
each other by a flop of one of the P^s, see figure |2j The appropriate geometry underlying 
our solution is the symmetric resolution. Indeed, when \Q\, \fj,\, \e\ < 1, the equation A3. 18|) 
describes the mirror of the symmetrically resolved orbifold, with Q, fi, e 2 being exponentials 
of flat coordinates of the Kahler moduli space, as we discuss in more detail below. 

We also note that on general grounds it is known that for Calabi-Yau manifolds of the 
form uv + H(x,y) = 0, with H(x,y) = encoding a Riemann surface as in ( 13.1 8ft . the 
special geometry relations reduce to 

where A is a reduction of the holomorphic three-form along u, v directions, a and b are dual 
one-cycles on a Riemann surface H (x, y) = 0, and F op is the topological string free energy. 
The same relations hold for the free energy Fq of matrix models, if T is identified with the 
't Hooft coupling |27| . Therefore the fact that the spectral curve in the case we consider 
agrees with the mirror curve of Y, ensures the agreement of derivatives of matrix model 
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and topological string free energies with respect to T, up to an integration constant which 
is a function of Q and \x (which are just parameters of the matrix potential). As the exact 
topological string partition function is a symmetric function of Q, fi and e 2 , this implies 
that this integration constant must restore this symmetry, and the resulting matrix model 
free energy 

F = Lis(Q) + Li 3 (/i) + Li 3 (e 2 ) + Li 3 (Q/xe 2 ) - Li 3 (Qe 2 ) - Li 3 (/ie 2 ) - Li 3 (Q/i), 

has to agree with the topological string result F op . 




Figure 2: Two resolutions of the C 3 /Z 2 x Z 2 geometry, symmetric one (a.k.a. 
closed topological vertex, left) and asymmetric one (right), related by a flop of 
one P 1 . The Kahler parameters of both geometries are related to each other 
as P x = Q 1 Q 2 ,P 2 = l/Q 2 ,P 3 = Q 2 Q 3 - 



For the BPS counting problem, we are interested in the limit of T — > oo, or equivalently 
e — > 0. With appropriate shifts of x and y, the equation (13. 18[) in this limit becomes 



li e 2y + e x+y + e x + (l + Qfi)e v + Q = 0. 



(3.20) 



The manifold Y corresponding to this curve is the SPP geometry, with Q and [i being 
exponentials of flat coordinates representing sizes of its two P^s, which encode two copies 
of the initial 0(—l, —1) — » P 1 geometry, see figure [21 Not only does the spectral curve 
agree with the mirror curve of the SPP geometry in the limit of g s — > 0, but in fact the 
matrix integral reproduces the full topological string partition function at finite g s . Indeed, 
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it is known that the SPP topological string partition function, with Kahler parameters Q 

and /i, is equal to 

[l - Qq k ) k {l - fiq k ) k 



z% p (i,q, ») = ![■£ 



(3.21] 



fc=i 



q k ) 3k / 2 (l - {iQq k ) k ~ 

On the other hand, from the explicit structure of the BPS generating function and formulas 
f !2.10p . ( I2.12p and (I2.13J) . we find that the value of the matrix integral, in the N — > oo limit, 
is related to the above topological string partition function as 



Z ma trix(g, Q; n) = Zf^(q, Q,fi = Q~ l q n ) ■ \J(l ~ Q 



k\k/2 



(3.22) 



fe=l 



In this way, the Kahler moduli Q and the chamber number n for the BPS counting on 
the conifold are unified into the two Kahler moduli of the SPP geometry. We note that 
there is an extra factor of the MacMahon function in this relation. The appearance of the 
MacMahon factor, which is independent of Q, is a common and subtle issue in relations 
between the topological string and other systems. 





Figure 3: Toric diagram for the Suspended Pinch Point (SPP) geometry, 
and the corresponding dual diagram. This manifold contains two copies of 
0{— 1, —1) — > P 1 geometry. 



We note that the spectral curves (13. 18[) and (I3.20J) arising from the matrix model 
automatically encode the relevant mirror map. For example, in the parametrization of 
( I3.20p . Q and \x are directly identified with the exponentials of the flat coordinates. We 
can verify this by explicit evaluation of period integrals, see [29], section 3.3. The form 
(I3.20p of the curve factorizes for jjl — 1 and Q — 1, which is consistent with degeneration of 
the topological string partition function A3 . 2 1 [) for these values. Also in the limit x — > ±oo, 
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the solutions of the curve equation for e y reproduce appropriate locations of the asymptotic 
legs of the SPP toric diagram in figure El The same parametrization naturally arises also 
in [30] as a characteristic polynomial of the dimer model; see the example in section 4.2 
and in particular (4.2.8) of [31J. All these arguments can be extended to the Y mirror 
curve (13.1 8p . Note that the standard parametrization of this mirror curve, such as the one 
in [32], would suggest the equation e x+y + e x + e y + e 2 e 2x + fie 2y + Q = 0. This is however 
valid for large values of Kahler parameters, and consistent with (j3.18p . as in this regime 
the quadratic terms in Q, fj, and e 2 are negligible. 

Because of the form of the spectral curve at finite 't Hooft coupling (I3.18p . it is natural 
to conjecture that the partition function of our conifold matrix model for finite 't Hooft 
coupling is equal to the topological string partition function of the resolution of C 3 /Z2 x Z2 
[28], modulo a MacMahon factor 

-total , ni * = f^ n _k )k TT (1 ~ Q? fc ) fc (l ~ /*<ffll ~ <V)*(1 ~ Q^q k f 

(3.23) 

We chose the MacMahon factor in such a way that it reduces to our result (I3.22p in 
the infinite 't Hooft coupling limit e — > 0. As another evidence for the conjecture, we 
point out that, in the limit Q,fi — > 0, our model reduces to the Chern-Simons matrix 
model (discussed in the next section) and the above partition function correctly reduces 
to the appropriate Chern-Simons partition function. It would be interesting to test this 
conjecture, for example by applying matrix model recursion relations of |15] . 

As discussed in [28j, the right-hand side of ( I3.23P is precisely (including the correct 
power of MacMahon function) the generating function of plane partitions in a finite K x 
L x M cube, and up to one power of MacMahon reproduces the closed topological vertex 
partition function with Kahler parameters identified as Qi = g s K, Q2 = g s L, Q3 = g s M 
(this generalizes C 3 model of plane partitions with one wall discussed in section |27L|) . In 
the present case we have one analogous identification of the 't Hooft parameter T = g s N. 
Our 't Hooft parameter has also a nice combinatorial interpretation: finite A^ corresponds 
to matrices with N eigenvalues, which in the construction of our matrix models arise from 
truncation of products in (I4.18P to N operators T'. This translates to truncation of Young 
diagrams, which arise from slicing of the crystal model pyramid, to at most N rows, which 
is equivalent to considering a wall at location N. Therefore our present model involves one 
wall associated to finite 't Hooft coupling, the second parameter \i which involves finite n 
(which also measures a size of the crystal), and the third parameter Q which appears in 
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the matrix model potential in the same way as /i, however does not have a clear crystal 
interpretation. The cube model of |28j involves three symmetric walls and has the same 
generating function (I3.23p . It would be interesting to understand the relations between 
these two models in more detail. 

As the final remark, we note that there are three limits in which our full matrix model 
reproduces both the mirror curve, as well as the topological string partition function of the 
conifold. The first such limit //, Q — >■ brings us to the Chern-Simons matrix model and 
will be discussed in the next section. The second limit fi, e — > is just the commutative 
limit of the model with matrices of infinite size. In both these limits it is not surprising that 
the size of the conifold is identified respectively with 't Hooft coupling e~ T or the original 
Kahler parameter Q. However in the third limit Q, e — > we obtain the conifold of the size 
fi = Q~ 1 e~ T , which in fact means the Q vanishes however the chamber parameter r — > oo. 
It also corresponds to the commutative limit, and shows that for vanishing Q the role of 
the conifold Kahler parameter is attained by fi. This is in agreement with the picturesque 
identification of the conifold size with the length of the top row of the pyramid in the 
crystal melting model, and puts this identification on the firmer footing. 

3.3 Relation to the Chern-Simons Matrix Model 

We now discuss the commutative chamber n — > oo of the conifold model presented 
above. We show that it leads to the matrix model which is equivalent to the Chern- 
Simons matrix model, and these two models can be unified in a geometric way. By the 
Chern-Simons matrix model [T2l fl~3| |2~T] we understand the unitary matrix model with the 
Gaussian potential, as in ( 13. ip . and finite 't Hooft coupling T = Ng s . Including the shift 
( lA.ip arising from the measure, we write its potential as 

V cs = T\ogU--(\ogU)\ duVcs= ~ v g , (3.24) 

Our present model is also unitary and in the commutative chamber the derivative of 
its potential (I3.6P reduces to 

T - \og(U + Q) 
OuVn^oo = — . (3.25) 

We recall that our matrix model arises from rewriting the BPS generating function, which 
in the n-th chamber takes form (I2.10p . In the commutative chamber n — > oo the term 
M{Q~ V ) = rijfc(l — Q~ l q k )~ k is removed from that expression. On the other hand, in this 
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limit the prefactor (12.131) reduces to a single MacMahon function M(l) = Ylki^ ~ 1 k ) k ■ 
Therefore in the commutative chamber we find 



\\dui \\iiij -u k f ] J 



^' " i j<k k 



C 9s -'™( n !)_ 



The ratio on the left hand side is precisely the partition function of the Chern-Simons 
theory on S 3 , which is also reproduced by the Chern-Simons matrix model (I3.24p . The 
spectral curve of that model has genus zero and is identified with P 1 which arises from the 
geometric transition of the S 3 . The size of this P 1 is given by the (finite) 't Hooft coupling 
T. Now we find the model whose partition function is given by the same Chern-Simons 
partition function and its spectral curve has also genus zero, however our association of 
parameters is different. Instead of finite 't Hooft coupling parameterizing the size of P 1 , in 
our model 't Hooft coupling is infinite, while the size of P 1 is encoded in a fixed parameter 
Q deforming the unitary Gaussian potential as in (I3.25p . As an immediate check we notice 
that for Q = our potential ( I3.25P indeed reduces to (I3.24p . and for infinite T it reproduces 
Gaussian result for plane partitions (13. ip . The dependence of the potential Kwoo on the 
parameter Q is shown in figure HI As we approach the conifold singularity at Q — 1, it is 
interesting to observe the flattening of the matrix potential. It is known that the conifold 
singularity has to do with the flattening of the Coulomb branch moduli space |3"3"j I3"4"] . This 
indicates a connection between the matrix variable and the Coulomb branch variables. 

With both finite 't Hooft coupling T and finite Q, we find a unifying geometric view- 
point, again in terms of the SPP geometry, however now with Kahler parameters Q and 
e~ T . In this topological string limit the equations (I3.12p and ( I3.13P take form 



y/a+ + Q + ^a_ + Q = 2e~ T '\ (3.26) 

xA-K + Q) + \A+(a- + Q) = (V^+Vb)e T/2 , (3.27) 
and their solution is given by 

o± = -1 + (2 - Q)e 2 ± 2te v /{l-Q){l-e 2 ), (3.28) 

which leads to the following form of the resolvent 

, x 1, ( u+l + Qe~ T - y/(u + 1 + Qe- T f - A{u + Q)e 2 \ 
*(«)^ = ^ ^ ( 2e . T(u ~ Q) j • (3.29) 

The spectral curve which arises from this resolvent is again mirror curve of the SPP ge- 
ometry and reads 

x + u + xu + x i t JL_ + t £_ = . (3 .30) 

17 



■V n ^oo ((f) 




<P 






Figure 4: Matrix potential (without L4 . 1\) shift), —V n -*ao{ ( p) '■ 
Li 2 (— Qe~ v ) in terms of a variable u = e v . The blue plot represents the Gaus- 
sian potential with Q = 0. Increasing Q flattens the potential (red and green). 
At the conifold singularity, corresponding to Q = 1, the potential becomes flat 
(black plot). 



It is clear that e — > limit leads to the conifold geometry with the conifold of size Q. 
Finally, for Q = the resolvent j5jZgJ) 



1, /u + l-J(u + l) 2 - Aue~ T 
WW ^» = ^T l0g { 2u^ 

agreed with the resolvent of the Chern-Simons matrix model found in [T5j 127] , and the 
spectral curve reproduces the conifold mirror curve of the size given by the 't Hooft coupling 

x + u + xu + x 2 e~ T = 0. 



3.4 C 3 /Z 2 

A similar analysis as for the conifold can be performed for C 3 /Z2 geometry, for arbitrary 
chamber. Even though we do not repeat a matrix model derivation of the spectral curve 
for this case, we note that the relation to topological string theory is also immediate, and 
the relevant geometry Y for this case is the resolution of C 3 /Z 3 singularity shown in figure 



2 Instead of introducing TlogU term to the potential (|3.24p to get the standard Vandermonde deter- 
minant, the solution in [2TJ involves completing the square, which leads to a redefinition Where = VW!\ e ■ 
Due to a different sign of g s we also need to identify 't Hooft couplings as There = — t 27 • Taking this into 
account, our cut endpoints (|3.28[) with Q = also agree with those in [27] . 
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This geometry contains two P-'-'s of O(0, —2) type, and denoting its Kahler parameters 
by Q and /i, its topological string partition function reads 



Z^(q,Q,fi) = H(l - q k y^ 2 {l - Qq k )-\l - M k y k (l - ^Qq 
fc=i 

Therefore, from (J2.18P and (I2.19P we find in this case 

oo 

^ ma tnx(g, Q\ n) = zSf 3 (q, Q,H = Q-\ n ) ■ 11(1 - q k ) k ' 2 



k\-k 



(3.31) 



(3.32) 



fc=i 



This shows that the matrix model partition function (12.181) in the n-th chamber is equal 
to the topological string partition function for C 3 /Z3 with its two Kahler moduli given by 
Q and jjl = Q~ 1 q n , up to the MacMahon function as in (13. 22 p . This unified geometry Y 
contains two copies of the initial C 3 /Z 2 resolution. It would be interesting to check what 
geometry Y would arise for finite 't Hooft coupling, and whether it is consistent with the 
total matrix model partition function. 





Figure 5: Toric diagram for the resolution ofC 3 /Z*3 singularity, and the corre- 
sponding dual diagram. This geometry contains two copies o/C 3 /Z 2 resolution. 



3.5 General Toric Calabi-Yau Manifold 

Matrix model for a general toric manifold and in general chamber can be constructed 
in a similar manner, following fermionic approach of |6l [7], and then analyzed along the 
lines above. We do not present a construction in a general chamber which is technically 
much more involved, however we found explicit expressions for matrix models for a general 
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manifold X in the non-commutative chamber. These models are presented in section 12.41 
Nonetheless, we postulate that for arbitrary chamber, in the infinite 't Hooft limit, we 
should also find a toric manifold Y which contains two copies of X. For finite 't Hooft 
limit matrix model spectral curve would encode yet more general manifold Y. 

4 Derivations of the Matrix Models 

In this section we give two derivations of our matrix models. One derivation (sec- 
tion 14. ip uses the free fermion formalism, while the other (section 14. 2p uses a set of non- 
intersecting paths] 3 ] 

Both derivations are based on the following observation. Let us begin with the crystal 
melting model of [5J. Given a configuration of a crystal, we can slice the crystal by a 
sequence of parallel planes. On each slice, we have a Young diagram. The Young diagrams 
evolve according to the interlacing conditions, which is equivalent to the melting rules of 

m 

For C 3 , we have jl] 

. . . -< A(-2) -< A(-l) -< A(0) y A(l) y A(2) y ... (4.1) 

where we write A >~ \x (equivalently \x -< A) for two partitions A = (Aj) and fi = (/^j), if 

Aj = p«i + 1 or Aj = fii for each i. (4-2) 

For conifold in the non-commutative chamber [35], 

...-< A(-2) i A(-l) ■< A(0) y A(l) y A(2) y ... (4.3) 

where we write X y ft for A = (Aj), fi = (/ij) if A* y (jf, i.e. 

...>X l >fi i > Aj_i > /Ui_i . . . (4.4) 

We can also discuss more general toric Calabi-Yau 3-folds X without compact 4-cycles 
(see figure [H]). The (p, q)-web for X has \ vertices, where \ is the Euler characteristics of 
X. To each vertex we associate a sign Si = ±1 so that, 

SiSi+l — Si, (4-5) 



5 We have been informed by Mina Aganagic that there is yet another derivation of the matrix models 
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where the sign factor Sj = ±1 is defined in section 2.4. This means that if the local 
neighborhood of z-th P 1 represented by an interval between vertices % and i+ 1 is 0(—2, 0), 
then Si+i = Sf, if this neighborhood is of 0{— 1, — 1) type, then SV+i = —Si. There is 
a binary choice of overall signs Si', the type of the first vertex could be chosen as either 
Si — +1 or Si = — 1. This choice corresponds to the exchange of rows and columns of 
Young diagrams. Each choice gives rise to a matrix model potential, and they are related 
to each other by analytic continuation. 




Figure 6: Toric diagram for Calabi-Yau manifold without compact 4-cycles 
arises from a triangulation of a strip. There are N independent P 1 's with Kah- 
ler parameters Qi = e~ li , and x vertices to which we associate © and signs Si. 
Intervals which connect vertices with opposite signs represent 0( — l, —1) — > P 1 
local neighborhoods. Intervals which connect vertices with the same signs rep- 
resent 0(—2, 0) — > P 1 local neighborhoods. 



Given such Si, the interlacing conditions in the noncommutative chamber are given by 
. . . % 2 A(-2) 8 J A(-l) 5 A(0) y A(l) y A(2) ? . . . (4.6) 

where >-=>- and we extended the definition of Si (i = 1, . . . , x) to Si (i G Z) by periodic 
identification: Si +X = 5j. More general expression, applicable to any chamber, is given in 

[a Eg. 
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4.1 Derivation (I): Free Fermions 
4.1.1 Wall Crossing and Free Fermions 

The first derivation is based on the free fermion formalism developed in [6j [7] (see also 
[36J), which we now review briefly. 

The basic idea is as follows. We have seen that states in the crystal melting model 
are represented by Young diagrams. Since Young diagrams are represented by states of 
free fermion systems and evolutions of Young diagrams by vertex operators, the partition 
function is written as a correlator of fermions bilinears |U [53] . 

We give the resulting expression in the notation of [B]. For any toric geometry with- 
out compact 4-cycles, the generating function of DT invariants in the non-commutative 
chamber can be written as 

Z BPS (q,Q;n = 0) = (n + \n_), (4.7) 

where \£l±) are fermionic states which will be defined below. Moreover, we can introduce 
wall crossing operators W p to write expression in other chambersp where n p = m for the 
p-th P 1 and all other n's set equal to zero: 

Zbps (q,Q;n p = m, all other n = 0) = (fi+|(W p (l)) m |fi_). (4.8) 

In the remainder of this subsection we give explicit expressions for the states \Q±) and 
wall-crossing operators W p . 

We first define a vertex operator T±(x) at each vertex as 



■■ s '- +1 '^-r ± (x), rg=- 1 (^) = r , ±(x) 



[X 



where V±(x) and T' ± (x) are defined in appendix [B] These operators represent the evolution 
of Young diagrams \(t), and T + , T_, I~" + and T'_ are nothing but the evolution rules -<, >-, -< 
and >-, respectively [35J. 

Next, we consider a product of \ such operators T±(x) interlaced with \ operators Qi 
representing colors q t , for i — 0, 1, . . . , % — 1. Operators Qi, . . . , Q x -i are associated to P 1 
in the toric diagram (also defined in appendix [B]), and there is an additional Qq. We then 
introduce 

Zt (x) = If (x)Qxr£ 2 (x)Q 2 ■ ■ ■ rf*- 1 (a;)Q x _ir£ x (x)Q . (4.9) 



In this example, inserting wall crossing operators is equivalent to commuting vertex operators, which 
is proposed in [71136). 
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Commuting all Q^s to the left or right we also introduce 

A + (x) = (Q Q 1 ---Q x ^)- 1 Mx) = T^(xq)T^(^-)...T^( ^ ),(4.10) 

Qi QxQ.2 • • • 9x-i 

A_(x) = A4x)(Q Q 1 ---Q x _ 1 )- 1 = T s J(x)T s J(xq 1 )---T^(xq 1 q 2 q x ^ 1 ). (4.11) 
The states \Q±) are now defined as 

(fi+| = (0\...A + (l)A + (l)A + (l) = (0\...A + (q 2 )A4q)A + (l), (4.12) 

|n_> = A_(l)Z_(l)A_(l)...|0) = A_(l)A_(g)A_(g 2 )...|0). (4.13) 

It was shown in [6J that we have the relation (J4.7P under the following identification between 
qi parameters which enter a definition of \ft±) and topological string parameters Qi = e~ Ti 
and q = e~ 9s : 

qi = (SiS i+1 )Qi, q = q q 1 ■ ■ ■ g x _i. (4.14) 

In addition the wall-crossing operators are defined by 

W p (x) = (r^aOQir'iCaOQa • • -r'_ p (x)Q p ) (t^ +1 (x)Q p+1 ■ ■■T t ?- 1 (x)Q x - 1 r t f(x)Q ) 

(4.15) 
and the relation (J4.8P holds under the change of variables 

Q P = (S p S p+1 )q p q m , Qi = (SiS i+1 )qi for i ^ p, q = q q 1 ■ ■ ■ q x _ v (4.16) 

4.1.2 Matrix Models from Free Fermions 

Once the BPS partition function is written in the fermionic formalism, it can be turned 
into a matrix model upon inserting appropriately chosen identity operator in the correlator 
(USD: 

Zbps (q,Q;n p = m, all other n = 0) = (0 + | I (W p {l)) m \tt_). (4.17) 

The identity operator I is represented by the complete set of states \R)(R\ (representing 
two-dimensional partitions). Using orthogonality relations of U(oo) characters \Ri an d 
the fact that these characters are given in terms of Schur functions \R = s r{u) f° r u = 
(u\, u 2 , W3, . . .), we can write 

R P,R 



J ID D 



s P t(u)s R t(u)\P)(R\ 

P,R 



/ dc/ (n r -(^)i°>) (<°i n r +(^ _i )), (4.i8) 
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where dU denotes the unitary measure for U(oo), which can be written in terms of eigen- 
values Ui = e 1 ^ of U as, 

dU = Y[d<j) k JJ (e** - e^) (e~* - e"^) . 

fc i<j 

Having inserted the identity operator in this form into (I4.17P we can commute away T' ± 
operators and get rid of operator expressions. This leads to a matrix model with the unitary 
measure dU . In case of the non-commutative chamber all factors arising from commuting 
these r' ± operators depend on U\ and contribute just to the matrix model potentials. In 
other chambers additional factors arise which do not depend on Ui and therefore contribute 
to some overall factor C n (in a chamber labeled by n). 

Thus in general we write the DT generating function as a matrix model, up to the 
factor C n . In the non- commutative chamber, the integrand can be expressed in terms of 
the theta-product, 

oo 

®(U\q) = JJ(1 + Uq k ){l + f/-y +1 ), 

fc=0 

and in other chambers of certain modification thereof. 

We emphasize here that this fermionic method of constructing matrix models applies 
to any chamber for any toric Calabi-Yau 3-fold without compact 4-cycles. This includes, 
for example, chambers where the BPS partition function becomes a finite product [6J. 

One may ask if our construction of the matrix model is unique. One potential source 
of ambiguity is the location of the operator I. In (I4.17P we inserted the operator I on the 
left side of (Wi(l)) m . When inserted on the right, we find a seemingly different matrix 
model potential, for example, 



[1 + Uq k ) (i + [/-y +rMM 

+ QUq k - n )(l + Q- l U- l q k + n + l ) 



^^n „ six-.-.a . «■«» 



\k=0 

in the conifold with the same prefactor (I2.13p . This integral can be turned into 



""'/^(wn^). <^> 



by the simple change of integration variable, U -> q n+1 Q~ 1 U~ 1 . The resulting integral 
is similar to our original integral ( 12. lip , but the numerator and the denominator are 
exchanged. Then the matrix model (I4.20p can be derived by taking advantage of the 
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freedom of changing overall signs Si — > —S^s, which was mentioned in section 14.1.11 The 
two matrix models can then be related by analytic continuation in q, as described in the 
footnote 5 of [37J. Another possible ambiguity involves renaming U — > U^ 1 in ()4.18|) . which 
does not affect the form of the measure, however may affect the form of the potential. 

4.2 Derivation (II): Non-intersecting Paths 

In this section we give yet another derivation of our unitary matrix models, based on 
the non-intersecting paths. The fundamental observation here is that states in the crystal 
melting model, i.e., a sequence of Young diagrams satisfying interlacing conditions, can 
be equivalently expressed as a set of non-intersecting paths on an oriented graphjj Using 
the Linstrom-Gessel-Viennot (LGV) formula described in Appendix C, we can express the 
result as a determinant, which is an integral over a matrix whose eigenvalues are the height 
of the non-intersecting paths. This gives a multi-matrix model. We finally simplify the 
matrix model into a 1-matrix model. 

4.2.1 C 3 

Let us begin with the simplest example, C 3 . Let us define 

h k (t) = \ N - k+1 (t) + k - 1, (4.21) 

for k — 1, . . . , N. Since X(t) is a partition, we have 

h k (t) < h k+ i(t), (4.22) 

for all t. We also have the boundary condition, 

h k (t) = k-l when \t\ large. (4.23) 

Moreover, (14.1 j) means we have, for each step t, 

h k (t + 1) - h k (t) = or -1, 

for t > and 

h k (t + 1) - h k (t) = or 1, 



5 The matrix model for C 3 is constructed recently by [TT]. We will generalize those arguments to conifold 
later. Also, even for C 3 the explicit expression of the potential for the 1-matrix model (|4.35|) in our paper 
seems to be new. 
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for t < 0. For later purpose, it is convenient to introduce the minus of the sign function 
a(t') (t'GZ+|), 



a(t')=l r (4.24) 




so that 

h k (t + 1) - h k (t) =0 or a(t + h, (4.25) 

Summarizing, we see that Young diagrams {A(£)} are equivalently expressed by a set of 
coordinates hk{t), which satisfy the conditions above. 

We can represent these coordinates as a set of non-intersecting paths on an oriented 
graph shown in figure [3j The coordinates of the k-th path at time t (specified by the 
t-th dotted line) is given by an integer hk(t). It is easy to see that there is a one-to- 
one correspondence between the coordinates h^it) satisfying the conditions above and the 
non-intersecting paths on the oriented graph. The inequality (I4.22p is translated into the 
non-intersecting condition. The step condition (I4.25P corresponds to the fact that we have 
two arrows for each vertex on the oriented graph, one with the same coordinate and another 
with coordinate increasing/decreasing by one unit. Thus, the BPS partition function can 
be expressed as a sum over non-intersecting paths, 



{hi(t)} : non-intersecting t 

where paths are assumed to satisfy the boundary condition (I4.23p . 
By the LGV formula (appendix [C]), (14.261) is equivalent to 



(4.26) 



■Xo) = / II dh ^ deti 'i ( G ^' i> *)) ' (4-27) 



(4.28) 



where the Green function G(i,j; t) is given by 

G(i,j;t) = q^ h <V 

x [5 (hi(t + 1) - hj(l)) + 6 (hi{t + 1) - hj(t) + a(t + 1/2))] . 

The discrete coordinates hi(t) are turned into continuous variables. The delta functions 
enforce the condition f !4.25p . The contributions including off-diagonal components of Green 



6 This oriented graph arises from the lozenge tiling of the plane, which is another way of representing 
crystal for C 3 . The paper [TT] uses this tiling to construct an oriented graph. For the derivation of matrix 
model in this paper, however, we do not need to invoke the notion of lozenge tilings. 
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t=-3 t=-2 t=-1 t=0 t=1 t=2 t=3 





t=-3 t=-2 t=-1 t=0 t=1 



Figure 7: Top: an oriented graph for C 3 . Middle: an example of 3 non- 
intersecting paths shown in red. The location of the k-th path at time t gives 
hk(t). Bottom: The corresponding evolution of Young diagrams. 
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functions correspond to intersecting paths, which cancel out by the sign of the determinant. 
We also need to set the boundary condition ( I4.23P 



hi(t)e{l,...,N}, |t|>l, (4.29) 

where N is an integer which we take to infinity at the end of the computation. Keeping 
N finite corresponds to taking only the first N paths. 

The delta functions can be generated by introducing Lagrange multipliers (pit') (£' G 



2/' 



-^ I ' d(j){t') e -T*vAW)) det ( e ^(*'+|)^(*')) det ( e -^(t'-|)^(0; 



det 



SMi! + \)- hjtf - \)) + 8{h{t' + \) - h 3 {t> - \) + a(t')) 



where the potentials Vt>(<t>{t')) depend on the signs of i! and are given by 

e -M<« = i + e '<MO (4.30) 

We can enforce the boundary condition (I4.29P by limiting — M < t < M for the height 
function h{t) and — M — | < t' < M + | for the Lagrange multiplier 4>{t') for sufficiently 
large M , and by introducing the factors, 

A (V#(-m-§)) A (>( A '+I)) , (4.31) 

at the initial and final points, where A(e J?i ) is the Vandermonde determinant, 

A(e^) = det« {e tk ^) . 

In the free fermion formalism in the previous subsection, the two determinants in f)4.3ip 
correspond to the bra and ket states for the Fock vacuum. We also set the potentials at 
the two end points to vanish, V t , = _ M _i = V t , =M+ i = 0. The partition function then takes 
the form[] 

Z matti M = f\[dh{t) !\{d<P{t') A (e~^- M -^) A (e^ M+ ^) , 

* *' ' ' (4.32) 

x ]J g TrMt)JJ e -T r y t , W 0) det ( e ^'+|Wt')) det ( e -^'-|WO). 



t 



7 We drop an overall constant here for simplicity. In the final expression, we will present the correct 
formula including the overall factor. 
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We can turn this into a matrix integral. We use the Itzykson-Zuber integral over unitary 
matrix U [351 135| 

det(e^) = A(X)A(F) f dU e iTrXUYU \ (4.33) 

to generate squares of the Vandermonde determinants A(h(t)) 2 and A(0(£')) 2 for all t and 
if, except for <j>(-M - ~) and <f>(M + §), for which A(0(-M - ±)) and A(0(M + |)) are 
generated. The resulting expression can now be written as a matrix integral 



<{q)=jY[dH(t)Jl[d*{lf) 



A e 



i^(_Af_i)\ A /^(M+|) 



./ -,;•' A(0(-M-|))A(0(MH-|)) (4.34) 

x TTg^ H(t) FT e_Tr yt ' (<J>(<,)) e Tr i$ (*')Wt'+2)--ff(*'-|)) ) 
t f 

where H{t) (t G Z) and $(£') (t' G Z+ |) are N x N Hermitian matrices whose eigenvalues 
are (hi{t)) i=1> _ )N and (<pi{lf))i=i,...,N, respectively. 

The matrix model (j4.34p is a multi-matrix model. However, H(t) appear only linearly 
in the integrand, and can be trivially integrated out, yielding the constraints 

$(£') = $(£' -1) + log g-l iVxW = $f--j + (t' + -J\ogq-l NxN . 
The matrix model then simplifies to a one-matrix. The measure factor for $ = $(— |) is 



rf$- 



Me**) 2 



dC/, 



A(0) 

where U = e**. The factor (j4.3ip we inserted to impose the boundary condition turns the 
hermitian measure for $ into the unitary measure for U. Thus, the matrix integral can be 
written as 

Z matlix = JdUQ(U\q). (4.35) 

This gives another derivation of (12. ip . 

4.2.2 Conifold and More General Geometries 

Let us next describe the conifold in the non-commutative chamber. Again, we define 
hk(t) by (14.2 ip . We then have the non-intersecting condition (I4.22p and the boundary 
condition (14.23)) . We also have, from ( 14 . 3 )) . 
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1. When t is odd, 



2. When t is even, 



for t > and 



for t < 0. 



h k {t + 1) - h k {t) =0 or a{t + -). (4.36) 



< h k -i{t + 1) < /*(*) < h k (t + 1) < /i fc+1 (£) < . . . . (4.37) 



< h k -!(t) < h k {t + 1) < h k {t) < h k+l {t + !)<.... (4.38) 



The set of coordinates h k (t) satisfying the conditions above can be expressed as coordinates 
of non-intersecting paths of the oriented graph shown in figure [HJ Let us show that this 
is indeed the case. When t is odd, the story is similar to the C 3 example; we have two 
possibilities in (I4.36p . which corresponds to the two arrows starting from a vertex of the 
oriented graph, one with the same height and another with increasing/decreasing height 
by one unit. The situation changes when t is even; during one time unit, after going one 
unit horizontally we can choose to go vertically as much as we want, as long as we respect 
the non-intersecting condition. In other words, the intersection of the k-th path and the 
time slice t is an oriented interval, and when take h k (t) to be the value at the endpoint 
of the oriented interval, the interval is expressed as [h k (t), h k (t + 1)] (this is for t < 0; for 
t > 0, we have [h k (t + 1), h k (t)] instead). This means that the conditions f)4.37p . (I4.38P are 
translated into the non-intersecting conditions for paths. In the vertex operator formalism 
explain in section (14.1. 2ft . t odd (t even) corresponds to T± (T' ± ). 

The procedure to obtain the matrix model is similar to the C 3 example. The multi- 
matrix model is given as follows: 

w 9 ) =/n^ w /n»(o A k ( _M-j))A(U + i» (439) 

x TT <h {t) W e _Tt Vt ' m ' )] e^ **(*') W*'+|)-ff(t'-i)) ) 
t t> 

The main differences from the C 3 example are that we have two parameters which depend 
on time t as 

J g (fcodd), 
q t = I (4.40) 

I gi (t even), 
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\A \A N \ 

AAA N f\ 

AAA JS JS 

^i V! N K 



t=-3 t=-2 t=-1 t=0 t=1 t=2 t=3 



I/: I/: N \ 



t=-3 t=-2 t=-1 t=0 t=1 t=2 t=3 




t=-3 t=-2 t=-1 t=0 t=1 t=2 t=3 



Figure 8: Top: an oriented graph for the conifold. Middle: an example of 3 
non-intersecting paths on the graph shown in red. Bottom: the corresponding 
evolution of Young diagrams. 
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and that the potential takes the form 

e -VA4>) = ) K J, (4 41) 

[ 1/(1 -e^(*')) (£':even). 

When £ is odd, there are 2 possibilities for hk(t + 1) — hk(t), which is the reason for the 2 
terms in the potential. When t is even, e~ v *'^ = 1 + e - ** "^ ' + e -2 ^""^ ) + ...; this reflects 
the condition (which is part of (14.3 



h k {t + l) = h k (t)+ ma(t + 1/2), m = l,2,3 



i.e., hk{t + 1) > hk(t). The remaining conditions of (I4.38P is taken care of by the non- 
intersecting condition. 

Again, we can integrate out M t 's, and the matrix model simplifies. When we diagonalize 
it, we finally have 

W = /dCTdet|g|, (4.42) 

where U = e l ®^~i\ q = go<?i and Q = —q\. The sign in Q comes from the localization |40| . 
which should properly be taken into account in the definition of the matrix model. 

We can also repeat the same analysis for more general geometries. In particular, in the 
non-commutative chamber we obtain the results presented in section 12.41 The oriented 
graph can be constructed from the data of Si, by combining the 4 basic patterns (corre- 
ponding to T + , r_, r' + , r'_) for each «GZ; see figure |9j As an example, the oriented graph 
for SPP in the noncommutative chamber with S± = +1, 5*2 = — 1, 5*3 = —1 is given in figure 
IHJ We stress that this approach is equivalent to the fermionic picture described earlier. 
For example, the sum over all possible paths in the region t < (or t > 0) is encoded 
in the state (Q+\ (respectively |^-)). These states live in the Fock space associated to 
t = 0, and can be expressed in terms of a sum over two-dimensional partitions from both 
fermionic and non-intersecting paths viewpoints. The correlator ( 14. 7 p represents gluing 
paths extending in the t < region with paths in the t > region in a consistent way. 
Since the evolution rules of Young diagrams in more general chambers are already given 
in [61 [71 [36], it is in principle straightforward to generalize the analysis to more general 
chambers. 
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Figure 9: Top: An oriented graph in general cases are constructed by combining 
the 4 basic types of graphs shown here. Bottom: An oriented graph for SPP, 
with Si = +1, S2 = —1, S3 = — 1. 
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5 Discussion 

In this paper, we derived unitary matrix models of infinite-size matrices, which give the 
counting of BPS bound states of DO and D2-branes bound to a single D6-brane wrapping 
a toric Calabi-Yau manifold X without compact 4-cycle. These matrix models depend on 
a set of parameters Q, which keep track of the BPS charges, and the chamber parameters 
n. Both Q and n are associated to the Kahler moduli space Ai(X) of X. It turned out 
that these matrix models define the topological string on another Calabi-Yau manifold Y, 
whose moduli space contains two copies of Ai{X). The parameters Q and n are unified as 
the Kahler moduli of Y. In addition, when the 't Hooft coupling g s N is finite, we found 
yet more general manifold Y. In the crystal model this finite 't Hooft coupling has an 
interpretation of restricting a crystal configuration by a wall located at position N, and 
then the limit N — > oo provides mathematically rigorous definition of our models. 

The relation between the BPS counting on X and the topological string on Y is clearest 
in the commutative and the non- commutative chambers. In other chambers, there is a non- 
trivial prefactor in the relation between the BPS partition function and the matrix model 
partition function. We hope to understand the origin and the nature of the prefactor 
better. 

Our methods provide a rigorous derivation of matrix models and spectral curves, which 
encode the mirror map expected from the remodeling conjecture [14J. In this context it is 
interesting to note the subtlety related to the counting of MacMahon factors. For example, 
in the conifold example in the commutative chamber with either Q = or e _T = 0, we 
have one power of MacMahon function M(q), which agrees with topological string result 
and Chern-Simons partition function. However there is a mismatch by M(g) 1//2 between 
our matrix model integral formula and the topological string partition function for SPP. 
Similar mismatches arise in matrix models derived in \TT\ Q2J fIT\ . 

The notion of the spectral curve also exists in the dimer model. In |30], which discusses 
the thermodynamic limit of the crystal melting model, it was proven using the results 
of |41] , that genus contribution of the DT partition function in the noncommutative 
chamber agrees with the genus part of the topological string on the spectral curve of the 
dimer model, which is the mirror of X. An interesting problem is to understand how the 
spectral curve of the matrix model is related to that of the dimer model. 

The holomorphic anomaly equations of topological string amplitudes can be interpreted 
as the manifestation of their background independence [4"2| 145] . The relation between 
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the BPS partition function on X and the topological string on Y suggests that the wall 
crossing phenomenon on X may be related to the background independence on Y . In this 
context it would also be interesting to relate our analysis directly to the continuous limit 
of Kontsevich-Soibelman equations j44j . 

In this paper we considered bound states of D6-D2-D0 branes. This analysis can be 
extended, both from M-theory and matrix model, to include an additional D4-brane and 
associated open BPS invariants |19| 150] . Refined versions of our results can also be found 
using similar techniques [51| . 
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A Unitary Measure and the Migdal Integral 

Matrix models derived in this paper, either from fermionic or non-intersecting paths 
viewpoint, are of the form 



J matrix 



fdUe^ TTV ^^ u \ 



where the unitary measure, after diagonalization U = diag(tti, ...,u n ) with eigenvalues 
Ui = e"^, takes form 

dU = \\d<p k JJ(e* - e^)(e-* - e'^). 

k i<j 

This measure can be turned into the form involving the standard Vandermonde determi- 
nant dU — ¥ Ylk d u k rL<j( M * — u j) 2 a t the expense of introducing an additional term T log U 
to the matrix potential 

V unitavy (U) -)• V(U) = Knitary^) + T log U, T = g s N. (A.l) 

To find the resolvent for compact domain of eigenvalue distribution, arising from the initial 
unitary matrix ensemble, one can use results of |52j . Namely, the resolvent u(u) of the 
resulting matrix model can be solved using the Migdal integral, as also explained in [27] 
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and confirmed in explicit computations e.g. in [TSIH5] . In case of the one-cut matrix model 
this integral takes form 



^)4/- ^^ ^S , (A.2) 

v; 2TJ2iri u-z y/{z-a+){z-a-) 

where the integration contour encircles counter-clockwise two endpoints of the cut a±. 

In computing such Migdal integrals we often come across the situation where the deriva- 
tive of the potential d z V(z) contains terms of the form -2Ii£±£i_ j n this case we find 



1 j dz log{z + c) y/(u-a+)(u-a_) 

UJAU) = —= <t ; r- ; = = A. 3) 

cK J 2Tf 2m z{u - z) y/( z - a+)(z - a_) V J 

1 . / y / (a+ + c)(a- -u) - ^/(a_ + c)(a + - u) \ 2 
2uT V (u + c)(s/a_ — u — y/a + — u) J 

\/(u - a + )(u - a) ( a/(«+ + c)a_ - \/( a - + c ) c 
lor ' 



2«T A /a + a_ V c(y/a^~ — y/aZ) 

This result arises from contour integrals around poles at z = and z = u, as well as along 
the branch cut of the logarithm (—00, — c). To find the latter contributions the following 
integral is useful 



dx 1 (a/(x — a)(b — u) — \/{x — b)(a — u)) 2 

- log 



(x — u)\J(x — a)(x — b) \/( u — a)(u — b) (u — x)\J{u — a){u — b) 

In particular, for the conifold matrix model with the potential given in ( 13. 6p . the 
resolvent can be expressed as 

w(u) = u Qe -r{u) - uq{u) + — 1 —= + -). (A.4) 

21 V Uy/alaZ uJ 

In consequence we find that the resolvent is given by a sum of two terms, which in the limit 

u — > 00 are respectively constant and of order 1/u. Imposing the asymptotic condition on 

the resolvent u(u) ~ 1/u given in (13. lip implies that the constant term must vanish, while 

~ 1/u term must have a proper coefficient. This leads to the result (I3.10p . and moreover 

gives rise to the two equations (I3.12p and (13. 13[) for the endpoints of the cut a±. The 

solution to these equations is given in (I3.14J) . For various computations concerning this 

conifold example it is advantageous to use the identifies 

1-QeV 

1-^ie 2 ) ' 
l-Q(l-e 2 + /ie 2 )x2 



a + a - = 1 2 

v 1 — \ie l 



(o+ + Q)(o_ + Q) 

[1 + a + /i)(l + a_//) 



l-/ie 2 
1 - fi(l - e 2 + Qe 2 ) 



• 



1 — fie 



2 
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B Free Fermion Formalism 

For completeness we review free fermion formalism [46] following conventions of [6, 35J. 
We start with the Heisenberg algebra 

and define 

V±(x) = e^"> ^ a±n , r±(x) = e^">° « a±n . 

They act on fermionic states |/x) corresponding to partitions /i as 

r_(x)|/i> = l^xW-wiA), r + (x)|/i> = 5> H - |A| I*>. (B.i) 

A>-/x /x>-A 

ri.(*)l/*> = $> |AhH |A>, r' + (x)|/.) = J^x^-WIA), (B.2) 

A^/i /x>-A 

where >~ and >- are interlacing relations defined in (14. 2 p and ( 14. 4p . These operators satisfy 
commutation relations 

T + (x)T4y) = ^_r_(y)r + (x), (B.3) 

1 — xy 

r' + (x)r'_(y) = -J— r'_(y)r' + (x), (b.4) 

1 — xy 
T' + (x)T.(y) = (l + xy)T4y)r + (x), (B.5) 

r + (x)r'_(y) = (l+xy)r'_ (y)r + (z). (B.6) 

We also introduce various colors g 9 and the corresponding operators Q 9 

Q 3 |A) = gl A '|A>. 

They commute with T operators as 

T + (x)Q g = Q 9 T + (xq g ), r' + (x)Q g = Q g T' + (xq 9 ), (B.7) 

Q g T_(x) = T4xq g )Q g , Q s r'_(x) = V'_{xq g )Q g . (B.8) 

C LGV Formula 

In this appendix we explain the Linstrom-Gessel-Viennot (LGV) formula |4"7[l4"8]. which 
is crucial for the derivation of the matrix model in section 14.21 
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Consider an oriented graph without closed loops. We assume that a weight w(e) is 
assigned to each edge e of the graph. We consider N particles which follow paths pi, each 
starting at vertices Oj and ending at hi (i — 1, . . . , N). For such paths P = {pi : ai — > 6j}, 
we assign a weight 

w( Pt ) = l[w(e). (C.l) 

What we want to compute is the quantity 

F({a t },{b l })= Y, IJ^)' ( C - 2 ) 

P: non-intersecting i 

where the summation is over non-intersecting paths. LGV formula states that this can 
be computed by summing over general (meaning, including intersecting) paths. More 
precisely, when we define the "Green function" 

G(a t ,b,)= Yl W (P^ ( C - 3 ) 

p: a path from Oj to bj 

then LGV formula states that 

F({o i },{6 i }) = det(G(o i ,6 J -))- (C-4) 

The proof is elementary, and proceeds by checking that contributions from intersecting 
paths cancel out due to the sign in the definition of the determinant. The determinant in 
the formula can be thought of as a discretized version of a Vandermonde determinant for 
free fermions, representing the Coulomb repulsions among particles. 

Now consider a more general situation. Suppose that we are given a set of vertices 
{cti(k)}, where k = 1, . . . ,L. We consider N particles, with the following condition: 2-th 
particle starts from Oj(0), goes through Oj(l), then Oj(2), . . . , and finally arrives at Oj(L). 
Then the multiplicative property of the determinant says that 

L-l 

det(G(oi(l), aAL))) = TT det(G(a i (ib), a,- (A; + 1))) (C.5) 

fc=l 

This is the expression we need in the main text. 
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